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Abstract. 

Based on a conjecture about the diffusion tensor of rela- 
tivistic particles perpendicular to the magnetic field at a shock, 
and considering particle drifts, we develop a theory to account 
for the cosmic ray spectrum between 1 GeV and 10 4 GeV. The 
essential assumption is that the mean free path perpendicu- 
lar to the magnetic field is independent of energy and has the 
scale of the thickness of the shocked layer. We then use the 
basic concept that the cosmic ray particles are accelerated in a 
supernova shock that travels through the interstellar medium. 
Radio observations, with additional supporting evidence from 
optical and X-ray data, suggest that small additional and tran- 
sient bowshocks in locally free expansion contribute to the in- 
jection, while it is the overall shock which accounts for injection 
as well as the main acceleration of energetic particles, and we 
obtain a reasonable picture of the physical process with such a 
concept. Physically important ingredients besides the presence 
of a strong shock are diffusion, drifts, convection, adiabatic 
cooling, the injection history, and the topology of the magnetic 
field, here assumed for simplicity to be homogeneous in the in- 
terstellar medium. The result is a spectrum, which for strong 
shocks in a gas with adiabatic index 5/3 yields a spectrum of 
E~ 2A ' 2 . Interstellar turbulence with a Kolmogorov spectrum 
then leads by leakage from the galactic disk to a spectrum 
which is E- 2 - 75 , as observed in cosmic rays, and as deduced 
from radio observations of the nonthermal emission from our 
Galaxy as well as that of all other well-observed galaxies. We 
argue that the ratio of cosmic ray electrons to protons is deter- 
mined by the amount of expansion which takes place from the 
cessation of electron injection to the break-up of the shell by 
cooling instabilities. Since the highest particle energy reached 
derives from geometrical arguments, it depends on the charge 
of the nucleus and so higher Z elements are predicted to reach 
higher energies. 
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- supernovae - supernova remnants - shock waves 
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1. Introduction 

The origin of cosmic rays is still not completely understood. 
There are a few well-accepted arguments: 

a) The cosmic rays below about 3 10 4 GeV are due to the 
explosions of stars into the normal interstellar medium (Lagage 
& Cesarsky 1983); 

b) The cosmic rays from near 10 3 GeV up to the knee, at 

5 10 6 GeV, are very likely due to explosions of massive stars 
into their former stellar winds (Volk & Biermann 1988). The 
consequences of this concept have been checked by calculating 
the cosmic ray abundances and comparing them with observa- 
tions (Silberberg et al. 1990); the comparison suggests that up 
to the highest energy where abundances are known, this con- 
cept successfully explains the data. It is especially interesting 
that no mixing from the supernova ejecta is required to ac- 
count for the known cosmic ray abundances, Wolf Rayet winds 
and other strong stellar winds around evolved stars as sources 
are all that is needed at present; 

c) For the energies beyond the knee there is no consensus; 
Jokipii & Morfill (1987) argue that a galactic wind termination 
shock might be able to provide these particles, while Protheroe 

6 Szabo (1992) propose an extragalactic origin, although in ei- 
ther case the matching of the flux at the knee from two different 
source populations remains somewhat problematic; 

d) For the cosmic rays beyond the ankle, at about 3 10 9 GeV 
an extragalactic origin is required because of the extremely 
large gyroradii of such particles. 

In the first paper of this series (Biermann 1993, hereafter 
paper CR I) we have proposed that a major effect to consider is 
the transport coefficient in a geometry where the magnetic field 
is perpendicular to the shock normal, commonly referred to as a 
perpendicular shock, and where the shock is spherical. Then we 
investigate the propagation of a shock wave in either a stellar 
wind or into a homogeneous medium. In such a case particle 
drifts at the shock and in the upstream and downstream regions 
are important (see, e.g., Jokipii 1987). 

Observations can be a guide here: The explosions into the 
interstellar medium are explosions where the magnetic field is 
nearly perpendicular to the shock direction over most of 47r 
steradians. Radio polarization observations of supernova rem- 
nants clearly indicate what the typical local structure of these 
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shocked plasmas is. The observational evidence (Milne 1971, 
Downs & Thompson 1972, Reynolds & Gilmore 1986, Milne 
1987, Dickel et al. 1988b) has been summarized by Dickel et 
al. (1991) in the statement that all shell-type supernova rem- 
nants less than 1000 years old show dominant radial structure 
in their magnetic fields near their boundaries. There are several 
possible ways to explain this: Ray leigh- Taylor instabilities be- 
tween ejected and swept-up material can lead to locally radial 
differential motion and so produce a locally radial magnetic 
field (Gull 1973). It could be due to strong radial velocity gra- 
dients of various clumps of ejecta, or result from clouds being 
overrun that now evaporate and cool the surrounding material 
(all of which are mentioned by Chevalier [1977]). 

It may seem paradoxical that we will be considering per- 
pendicular shocks where the observations show a parallel mag- 
netic field, but it must be remembered that the observed field 
is a superposition of the configurations of all emitting regions 
in the telescope beam. The degree of polarization in all of the 
objects referred to above is low (typically 10% from a perfectly 
uniform field (generally 65% being both parallel and perpendic- 
ular components present. In young shells the former dominate, 
presumably reflecting the importance of the radial motions in- 
volved. (Although not relevant to the polarization, it might 
be noted that the parallel components will cancel on average, 
while the perpendicular components do not since they derive 
from the preexisting magnetic field structure.) 

The radial magnetic field configuration is also found in a 
few shell remnants somewhat older than 1000 yr, such as RCW 
86 (SN 185) and Pup A (with a kinematic age of some 3700 
yr, Winkler et al. 1988), showing that radial motion in some 
form persists beyond the youthful stage. In most "mature" 
remnants, however, the field is tangential, or complex (Milne 
1987). While a radially-directed magnetic field provides the 
most compelling evidence for large-scale differential motions 
parallel to the shock direction, there are other indications as 
well. It is generally believed that the oxygen-rich fast moving 
knots (fmk's) in Cas A are ejecta from the progenitor (Kam- 
per & Van den Bergh 1976). Oxygen- rich knots have also been 
seen in Pup A (Winkler et al. 1988) and the similarly aged 
G 292.0+1.8 (Braun et al. 1986), and may also be present in 
a number of extragalactic SNRs (Winkler & Kirshner 1985). 
Radially-directed radio structures have been mapped in Cas 
A, and their motion has also been determined (Braun et al. 
1987). Recent Rosat observations of the older (10 000 yr) Vela 
SNR show similar features in the X-ray emission which have 
apparently penetrated the main SNR shell (plenary talk by J. 
Triimper presented at the Texas/Pascos '92 Symposium, De- 
cember 1992) and are believed to result from fmk-like ejecta. 
(Their very soft X-ray spectra would make similar features dif- 
ficult to detect in more distant remnants.) 

The important conclusion for us here is that there appear 
to be strong radial differential motions in perpendicular shocks 
which provide the possibility that particles get convected par- 
allel to the shock direction. We assume this to be a diffusive 
process, and note that others have also pointed out that diffu- 
sion may be the key to understanding shock acceleration (e.g. 
Falle 1990). Our task here will be to derive a natural velocity 
and a natural length scale, which can be combined to yield a 
diffusion coefficient. A classical prescription is the method of 
Prandtl (1925): In Prandtl's argument an analogy to kinetic 
gas theory is used to derive a diffusion coefficient from a nat- 
ural scale and a natural velocity of the system. Despite many 



weaknesses of this generalization Prandtl's theory has held up 
remarkably well in many areas of science far beyond the origi- 
nal intent. We will use a similar prescription here. 

Consider the structure of a layer shocked by a supernova 
explosion into a homogeneous medium in the case that the 
adiabatic index of the gas is 5/3 and the shock is strong. Then 
there is an inherent scale in the system, namely the thickness 
of the shocked layer, in the spherical case for a strong shock 
limit r/12. There is also a natural velocity scale, namely the 
velocity difference of the flow with respect to the two sides of 
the shock. Both are the smallest dominant scale in velocity 
and in length; we will use the assumption that the smallest 
dominant scale is the relevant scale again to derive diffusion 
coefficients and other scalings. 

In paper CR I we introduced the basic theory on the basis 
of a transport coefficient as described above, and in paper CR 
II (Biermann & Cassinelli 1993) we tested it with the observa- 
tions of nonthermal radio emission from accelerated particles in 
shocks traversing stellar winds, while in paper CR IV (Stanev 
et al. 1993) we checked the spectrum and chemical composi- 
tion against air shower data in the particle energy range from 
10 4 GeV to 3 10 9 GeV. The interpretation of the radio observa- 
tions of stars before and after they explode led us to the concept 
that there is a critical shock velocity of order 10 9 cm s _1 above 
which electrons get injected into the acceleration process in a 
fashion similar to the nuclei, and below which this injection 
is severely restricted, an effect of several orders of magnitude. 
One possible choice for this critical velocity is that at which 
the post-shock thermal electrons reach relativistic velocities, 
or, in equipartition, the whistler velocity becomes relativistic. 

Biermann (1992), Rachen (1992), Rachen & Biermann 
(1992, 1993: paper UHE CR I) have proposed that the ultra- 
high energy particles arise from hot spots in nearby radio galax- 
ies; this hypothesis leads to a successful and nearly parameter- 
free explanation (Rachen 1992) of the intensity, chemical com- 
position and spectrum of these particles (see Rachen et al. 
1993, paper UHE CR II) with the important proviso that 
the mean free path in intergalactic space should be not much 
smaller than the characteristic distances between the sources 
and us, and may be similar to the scale of the large-scale bub- 
bles in the universe. These results independently confirm the 
existence of a steep spectrum component of heavy nuclei, iden- 
tified with galactic sources in paper CR IV. 

The observations of the cosmic rays themselves, and of the 
nonthermal radio emission from our Galaxy as well as from all 
other well-observed galaxies strongly suggest that in all care- 
fully studied galactic environments, the cosmic rays have a 
universal spectrum of very nearly E~ 8 ^ 3 below the knee at 
510 GeV (below about 10 GeV for the energetic electrons); 
direct air shower experiments show the spectrum beyond the 
knee to be well-approximated by E~ 3 . This overall spectrum 
is clearly influenced by propagation effects, since particles at 
different energies have different probabilities of escaping from 
the disk of the galaxy. It appears to be a reasonable hypoth- 
esis to approximate the interstellar turbulence spectrum by a 
Kolmogorov law, which leads to an interstellar diffusion coef- 
ficient proportional to E 1 ^ 3 . Such an energy dependence then 
requires a source spectrum of Cosmic Rays of approximately 
_E~ 7 / 3 below the knee, and approximately E~ 8 ^ 3 above the 
knee. 

In this paper we propose to derive such a spectrum for 
particle energies below about 10 4 GeV. The basic hypothesis 
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is, again, that we consider explosions into a homogeneous in- 
terstellar medium. We will use the standard Sedov solution to 
the expansion of such an explosion into the interstellar medium 
(e.g. Cox 1972). 

2. Definition of the task 

We consider the acceleration of particles in a shock that prop- 
agates into a homogeneous magnetic interstellar medium. The 
acceleration of particles is governed by the standard theory 
(Parker 1965), which includes the effects of first, time change, 
then radial diffusion, latitude diffusion, radial drift, latitude 
drift, then compression and finally sources (see paper CR I). 
The orientation which we use is such that the symmetry axis 
is parallel to the prevailing magnetic field, and we assume that 
to first order the magnetic field structure is unidirectional and 
homogeneous in the undisturbed interstellar medium. We take 
the cosine [i of the colatitude as our coordinate along 9. The 
interstellar magnetic field is known to contain both systematic 
ordered components and a random component (Beck 1986) of 
similar strength; hence we will demonstrate below (Sect. 8) 
that the simplification of using a uniform field does not in- 
troduce any unnecessary restrictions. The components of the 
diffusion tensor of interest here are the radial diffusion term 
K rr and the latitude diffusion term Kee- 

The Sedov similarity solution for an explosion into a homo- 
geneous medium can be written as a function of shock radius 
r and shock speed Ui; as a function of time £4, in units of 10 4 
years; explosion energy E51, in units of 10 5 ergs; and density 
of the interstellar medium n in units of cm" 3 . This gives (Cox 
1972) 



r = 13.7t 2 J 5 (E 51 /n ) 1/5 pc 
and 

Ui/c = 1.810- 3 t- 3/5 (S 5 i/n o ) 1/5 

= 0.090 (r/pc)- 3/2 (£ 5 i/n ) 1/2 . 



(1) 



Among the modifications which can be made to Eqs. (1) and 
(2) with practical application to SNRs, it is possible to account 
for (radial) density gr adien ts in both the ambient medium and 
the ejecta. Chevalier (1982) suggests, for example, that a mod- 
ification which may be appropriate for Tycho would result in 
r oc t 4 ^ 7 instead of the 2/5 power-l aw in Eq. (1). Kinematic in- 
vestigations of Tycho (Strom et al. 1982 ), Kepler (Dickel et al. 
1988a) and SN 1006 (Long et al. "~ 
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are deviations from the ideal 2/5 power- law, they are not sub- 
stantial, and it is moreover unknown what their precise cause 
is. We shall consequently adopt the Sedov solution for the sake 
of simplicity and because it seems to be a fair representation 
of reality. 

Our boundary conditions are the usual: We inject particles 
at some low energy which is assumed to be independent of all 
relevant properties of the problem, i.e. not dependent on radial 
distance, magnetic field strength or latitude. Downstream we 
assume the flow to take particles out of the system with the 
normal probability AU2/C (see, e.g., Drury 1983), where U2 is 
the downstream velocity relative to the shock. 

We propose to derive the essential properties of the parti- 
cle distribution function by analytical means, using heuristic 
arguments. The key will be the form of the diffusion tensor, 
especially the radial component n rr . 



3. Diffusion perpendicular to the average magnetic 
field 

Our basic conjecture (see CR I), then, is that the convec- 
tive random walk of energetic particles perpendicular to the 
magnetic field can be described by a diffusive process due to 
fast convective motions with a downstream diffusion coefficient 
K r r,2 which is given by the thickness of the shocked layer to- 
gether with the velocity difference across the shock, and is thus 
independent of energy: 



1 U 2 

fiVr.2 = 7T 77- r (Ui — U2) 

y Ui 



(3) 



Here we have assumed that in the combination of trans- 
port coefficients perpendicular and parallel to the magnetic 
field (see Sect. 8) the perpendicular convective transport dom- 
inates. As is obvious from Eqs. (51) and (52) below, this is 
a very good approximation for nearly all angles over the 4ir, 
since the parallel diffusion coefficient is linear with E, and the 
perpendicular convective transport is independent of particle 
energy. The exact limiting condition is (I — (j 2 ) / [i 2 3> Ej ' E max , 
which defines the very small solid angle over which the parallel 
diffusion is stronger. 

The upstream diffusion coefficient is then obtained (again, 
see paper CR I) by using the same velocity scale and the down- 
stream length scale multiplied by the density ratio, which gives 
r/3 for a strong shock. Again, the gyroradius of the same max- 
imum particle energy as in the shocked layer gives the same 
scale. Since the magnetic field is lower by a factor of U2/U1 
upstream, the gyroradius of the maximum particle energy up- 
stream is also r/3. Hence the upstream diffusion coefficient is 



«rr,l = g r (Ui — U2) 

It immediately follows that 

1 f 1 _Th 
9 V Ui 



Krr.l 



K-rr/2 



(2) rUi rU 2 



(4) 



(5) 



For these diffusion coefficients, it also follows that the resi- 
dence times (Drury 1983) on both sides of the shock are equal 
and are 



4 Krr.l 

Uic 



4 K rr ,i 

U 2 c 



4 r / _ Xh 
9c I Ui 



(6) 



Adiabatic losses then cannot limit the energy reached by 
any particle since they run directly with the acceleration time, 
both being independent of energy, and so the limiting size of 
the shocked layer limits the energy that can be reached to that 
where the gyroradius just equals the thickness of the shocked 
layer. This then leads to a potential maximum energy of 



£ m ax = l^ZerB 2 = \zerB 1 
6 Ui o 



(7) 



where Ze is the particle charge and Bi,2 is the magnetic 
field strength on the two sides of the shock. This means, once 
again, that the energy reached corresponds to the maximum 
gyroradius the system will allow, ft also means that we push 
the diffusive picture right up to its limit where the diffusive 
scale becomes equal to the mean free path and the gyroradius 
of the most energetic particles. We will derive a more restrictive 
condition below (Sect. 8). 
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4. Particle Drifts 

Consider particles which are either upstream of the shock, or 
downstream; as long as the gyrocenter is upstream we will 
consider the particle to be there, and similarly downstream. 

In general, the energy gain of the particles will be gov- 
erned primarily by their adiabatic motion in the electric and 
magnetic fields. The expression for the energy gain is well 
known and given in, e.g., Northrop (1963, equation 1.79), for 
an isotropic angular distribution where one term comes from 
the drifts and a second from the induced electric field (see pa- 
per CR I). We explicitly work in the shock frame, separate the 
two terms and consider the drift term first. The second term 
is accounted for further below, in Sect. 5. 

The #-drift can be understood as arising from the asym- 
metric component of the diffusion tensor, the #r-component. 
The natural scales there are the gyroradius and the speed of 
light, and so we note that for (Forman et al. 1974) 



K 8r — g fg C, 



(8) 



the exact limiting form derived from ensemble averaging, 
we obtain the drift velocity by taking the proper covariant 
divergence (Jokipii et al. 1977); this is not simply (spherical 
coordinates) the r-derivative of K,e r - The general drift velocity 
is given by (see, e.g., Jokipii 1987) 



E B 

C CUrle 777 • 



3Ze 



B 2 ' 



(9) 



The 6>-drift velocity is thus zero upstream and only finite 
downstream due to curvature: 



Vd,e = - — 
3 r 



(10) 



where r g is now taken to be positive (see above). 

It must be remembered that there is a lot of convective tur- 
bulence which increases the curvature: The characteristic scale 
of the turbulence is r/12 for strong shocks, and thus the cur- 
vature is 12/r maximum. Taking half the maximum as average 
we obtain then for the curvature a factor of 6/r which is six 
times the curvature without any turbulence; this increases the 
curvature term by a factor of six thus changing its contribution 
from 1/3 to 2 in the numerical factor in Eq. (10). Hence the 
total drift velocity is 



V d , e 



1 (lh\ erg 

2 \U 2 ) r ' 



(11) 



now written for arbitrary shock strength. It is easily verified 
that the factor in front of Eq. (11) is two for strong shocks 
where U1/U2 = 4. With AEi — we have then downstream 



AE 2 
E 



9 c V Ui) 



(12) 



which is the total energy gain from drifts. 

Let us consider then one full cycle of a particle remaining 
near the shock and cycling back and forth from upstream to 
downstream and back. The energy gain just due to the Lorentz 
transformations in one cycle can then be written as 



( E ) LT 



3 c 



(13) 



Adding the energy gain due to drifts, Eq. (12), we obtain 
AE 4(7i / U 2 



AE = Uh (,_U2\ 
E ~ 3 c V Ui) 



where 



x = 1 + 



(14) 



(15) 



5. Source expansion and particle injection history 

5.1. Expansion and injection history, shell 
Consider how long it takes a particle to reach a certain energy: 

-1 



dt 

Je 
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4Ui 
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Here we have used K rr ,i/U\ = K rr . 2 IU2. Since we have, 
from the Sedov solution [Eqs. (1) and (2)], r = (5/2) Uit, this 
leads to, 



dt 

T 



dE 3Ui 



2 K'pf \ 5 

Th x~rth 2 



E Ui 

and so to a dependence of 

E^ 



w = to (-) 



with 



P = 



3(7i 



Ui 



2 /vrr.l 5 

Th x rUi 2' 



(17) 



(18) 



(19) 



Particles that were injected some time t ago were in- 
jected at a different rate, say, proportional to r b . Also, in a 
d-dimensional space, particles have r d more space available to 
them than when they were injected, and so we have a correc- 
tion factor which is 



(!) 



i(b+d)f3 



(20) 



The combined effect [Eqs. (19) and (20)] is a spectral 
change by 



3Ui 2. K rr .l 



(21) 



We note that the factor 2/5 from the Sedov expansion drops 
out again. For a detailed argument on the sign conventions 
used here, see paper CR I. This expression can be compared 
with a limiting expansion derived by Drury (1983; eq. 3.58), 
who also allowed for a velocity field. Drury's expression agrees 
with the more generally derived expression given here. This 
also describes then the adiabatic losses, due to the general ex- 
pansion of the shock layer. Hence the total spectral difference, 
as compared with the plane parallel case, is given by 



3t/i 



Ui -U 2 



£(--1) + -(* + «0^ 

Ui \x J x rUi 



(22) 



If this expression is positive the spectral index is steeper. 

Furthermore, we have to discuss the effect of the expansion 
of the material in the shell, in this context already introduced 
by Drury (1983). This expansion is only that beyond the linear 
expansion of the overall flow already taken into account. We 
define the velocity for the flow field U(£r) as a function of 
radius £r interior to the shock radius r by 
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um = v(0€^i. 

We thus have, e.g., V(l) 
field can be approximated by 



v(0 



e + 1 

3£ 8 + 5' 



(23) 

V'(l) = 3/4. This velocity 
(24) 



E~ 



(30) 



As shown by Drury (1983), this velocity field introduces a 
steepening by 



3t/l Hrr,2 V'(l) 



Ui - U 2 rU 2 V{1) ' 



(25) 



neglecting drifts and in the limit of very small diffusion 
coefficient. However, in our context where we consider a finite 
shell, we have to use a non-local approximation. To do this 
we simply equate the expansion of the hot interior to the shell 
with the postshock speed of sound, which is given by 



C2 



27 



7- 



1 



1/2 



u 2 , 



(26) 



where the numerical factor is v5 = 2.23608 . . . for an adi- 
abatic constant of 7 = 5/3. The postshock gas expands at 
the shock in the observer's frame at Ui — U 2 , which is 3U 2 
for a strong shock. The linear expansion itself already provides 
(Ui — U 2 )/12, which is U 2 /4. Hence the extra expansion over 
linear is 2.75 U 2 — 2.23608 U 2 . Using then finite differences in- 
stead of derivatives, we have 



ni) 

V(l) 



2.751/2 - VEu 2 

3U 2 



gl2 
24 



1.970, 



(27) 



instead of unity as in Drury's case, which is 1.97007... 
times the term from the local flow field V(£). Here the factor 
23/24 comes from averaging the postshock radius from 1 to 
1 — 1/12, and thus describing the average of £/A£. Hence we 
have an extra term from this postshock adiabatic loss due to 
the shell expansion of 



3[/l 1 K rr . 2 V'(l) 



Ui-U 2 x rU 2 V(l) 



1.970, 



(28) 



which is then slightly less than double the effect which 
Drury discusses. For a strong shock then this particular effect 
adds 0.563 to the spectral index. 

For strong shocks in stellar winds, for which the natural 
length scale is r/4, the linear part of the expansion is three 
times higher and so the remaining difference to the postshock 
speed of sound is only 2.25 — VE, which gives an effect of only 
2 10 -3 in the spectral index, justifying a posteriori our neglect 
of this effect in papers CR I and CR II. 

The exponent b [Eq. (20)] describes the injection as a power 
of the radius; in a Sedov solution the injection, assumed to be 
proportional to pUi, is given by 



(29) 



Hence we have b — —3, so that b + d = 0. Thus we have for 
a pure Sedov solution, using the velocity field term introduced 
above, a spectral difference to the plane parallel case of 0.444 
and thus an injection spectrum of 



However, we note that this spectrum is valid for electrons 
obviously only if they are injected into the acceleration process 
at all, i.e. if the shock speed is above the critical velocity for 
electron injection. Adding the term for diffusive losses from the 
Galaxy the final observable spectrum expected is 



E~ 



(31) 



very close to what is observed, both in Galactic cosmic rays 
(protons and nuclei), as well as through the nonthermal radio 
emission of other galaxies (Golla 1989, electrons). 

Clearly, there is uncertainty in our method for treating the 
non-linear flow field; the various ways of correcting the term 
Eq. (25) all lead to numerical correction factors very close to 
2.0; above we have taken 1.97. Further non-linear corrections, 
such as also averaging the non-linear flow field itself, lead to 
spectral indices within 0.04 of the value given in Eq. (30). 

5.2. Expansion and injection history, fast moving knots 

Here we have to note that the radio images of Cas A (Anderson 
et al. 1991) suggest that there are many small additional bow- 
shocks which in fact may do some of the particle acceleration 
(see especially Braun et al. 1987). They are associated with 
compact clumps related to the fmk/ejecta (which, in fact, they 
may well be), and have travelled outwards essentially unde- 
celerated until their encounter with the high density material 
associated with the shock. Since the main shock has under- 
gone significant deceleration (Tuffs 1986, Braun et al. 1987), 
the clumps, expelled with a range of radial speeds, are now 
able to overtake it, brightening and hence becoming promi- 
nent as they do. Each clump is the apex of a bowshock which 
is radially oriented, and as it passes the shock front it becomes 
a fresh start into the interstellar medium and so can be pre- 
sumed to be in free expansion for as long as it survives, until 
it ultimately merges back into the overall general flow. The 
observations suggest that these bowshocks, initiated by small 
blobs of denser material, move with velocity constant in time, 
but different for different blobs; when such blobs happen to 
move faster than the present average velocity of the shell, they 
must have penetrated it in the past, but when their velocity 
is just the current value of r/t = (5/2)(7i, then they have just 
caught up with, and hence penetrate, the shell now. 

The drift velocity is then given by a consideration of free 
expansion into the homogeneous medium, and we approximate 
this free expansion here as a spherical shell of which only a seg- 
ment is realized. Writing for the blob radius r b and its velocity 
Uh we then obtain for the drift, again using the argument about 
convective scales [e.g. Eq. (11)], 

Vd ' B - 2UV 2 ~ b C - (32) 
The downstream residence time is then 
1 r b U h2 , TT T1 s 4 

Tr - 2 = n^^-^'fe' (33) 

The energy gain due to drifts is then given by 
E 9 c V U b J K ' 
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Now the argument above suggests that Uhi = (5/2) Ui. 
Using then the plausible assumption that ?7bi/fb2 = U1/U2 
we obtain 



AE = 5 Ui f 1 _Ui\ 
E ~ 9 c V U J ' 



(35) 



Adding this now to the Lorentz transformation energy gain 
(but also with the increased shock speed) we obtain for the 
parameter x again the value x = 7/6 [Eq. (15)]. In order to 
calculate the resulting spectrum we have to find now the ef- 
fective diffusion coefficient and express it in terms of the shell 
parameters, such as shell thickness of the overall explosion, for 
strong shocks r/12. For free expansion the parameters b and d 
are b + d = 3, and the expansion is linear so that the velocity 
field is V'(£) = 0. Hence we have to fix the residence time now 
in order to determine the spectrum. Now the total scale of the 
blob-induced convective turbulence is, upstream plus down- 
stream, (5/4)rb/3. This entire scale will grow until at some 
point it reaches the overall shell scale r/12; then by virtue of 
the switch-over in scales the overall shell scale will begin to 
dominate and the blob expansion will cease to be important. 
This means that 



rb = 



and also that 



r b ^bi = 2 rUl - 



(36) 



(37) 



Correspondingly we have to insert K rr ,i/(f(7i) = 1/24 and 
so obtain for the spectral difference from the plane parallel case 
0.715, which gives an injection spectrum of 

£- 2 ' 715 , (38) 

nearly 1 /3 steeper than the spectrum for the overall shock 
acceleration, Eq. (30), suggesting that these smaller shocks pro- 
vide just some injection for the larger shocks. However, the 
smaller shocks may not only run into a geometric limit, they 
are likely to run into an energetic limit even earlier (i.e. no 
longer have sufficient energy to be in free expansion locally), 
and so cease to be important before reaching the limit derived 
above. 

As can be seen from the expression for the spectral index, 
it depends on the blob radius, and thus on expansion time. The 
spectrum steepens with blob radius rb, which is proportional 
to expansion time for as long as the blob is in local free expan- 
sion. The spectrum implied is E~ 2114 at 0.3 r b , then E~ 2 ' 715 
at exactly rb, to steepen to E~ 4 ' 429 at 3rb. Since energetic 
particles experience both the effect of the small blobs as well 
as that of the overall shock, the overall spectral index for the 
blobs is expected to be slightly steeper than the average (non- 
blob) spectral index of the shell, consistent with observations. 
Wc thus surmise on energetic grounds that these blobs break 
up below the critical blob size given above, and note that their 
radio emission does in fact fade rapidly as they pass the shock 
front (Braun et al. 1987), which is consistent with the breakup 
requirement. This means that the blobs are likely to change 
over to a local Sedov expansion before reaching the critical 
blob size; this then would lead to a cessation of electron injec- 
tion as soon as the the shock velocity drops below the critical 
velocity. 



The spectrum derived is the one injected by small spherical 
blobs in local free expansion. This spectrum has more power 
proportional to the square of the shock velocity, since the shock 
velocity Uhi is 5/2 higher than the overall shell shock veloc- 
ity Ui, and also the effective area is, of course, approximately 
half a sphere, giving another factor of two over the area of a 
planar circle on the large shell. Hence the multitude of these 
little blobs puts about one order of magnitude more power 
per area of the large shell into the injection than the spherical 
shell shock directly. We conclude that the small blobs driven 
by the fast moving knots may contribute substantially to the 
injection of particles, just as surmised by Braun et al. (1987). 
(Parenthetically, it should be noted that the possible role of 
fmk's in particle acceleration was first investigated by Scott & 
Chevalier [1975] who, however, assumed a second-order Fermi 
process in the turbulent fmk wakes.) 

We add that these small bow shocks are not the only ob- 
jects that should have time-variable radio emission; radio ob- 
servations suggest that there are also many nearly stationary 
magnetic substructures, with likely a magnetic field up to a 
factor 5 above that given by the shock-enhanced interstellar 
field of about 20 fj,G. Such substructures can arise when the 
overall expansion overruns preexisting features in the inter- 
stellar medium, for which a straightforward analogy with the 
solar photosphere would suggest that magnetic substructures 
of much stronger magnetic fields might be quite common. 

There is one immediate check: Braun et al. (1987) and An- 
derson et al. (1991) emphasize that the large bow shock fea- 
tures with the steepest radio spectral index have the strongest 
time variation of their flux density. In our simple concept this 
can be interpreted as the consequence of two effects: First, with 
the size of the blob or bow shock the region emitting scales as 
the third power of the radius of the blob and so steadily in- 
creases. Second, with the steeper spectrum the weight of the 
particle spectrum is shifted to low particle energies and so the 
radio emission is increased. For going from a radio spectral in- 
dex of, e.g., —0.6 to —0.7 the combined effect is about a factor 
of 4, equally caused by both effects. 

The spectrum derived is very close to the observed galactic 
cosmic ray spectrum at low energies and we propose that this 
is the mechanism to generate it. 

6. The time evolution of the radio spectral index 

Since the shock velocity is steadily decreasing with time for a 
Sedov expansion, we have to consider the effect of particle ac- 
celeration through Alfven wave scattering, commonly referred 
to as second order Fermi acceleration. The importance of this 
mode of acceleration has been appreciated for some time (e.g., 
Scott & Chevalier 1975, Chevalier et al. 1976). 

The energy gain through this effect can be included in our 
derivation of AE/E: 

(¥).-(¥)„(£)'-- -(¥)„'• «> 



This translates into a time- dependent value of x: 
x = 1.167+ (r/r s f. 



(40) 



This in turn translates into a spectral index for the particle 
distribution which gets flatter with time and expansion radius: 
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1 + 



1.657 



1.167 + (r/r a ) 3 ' 



(41) 



with time if the efficiency of injecting electrons, rj, is a constant. 
The luminosity is then given by 



The radio spectral index of the nonthermal synchrotron 
emission then has a radius-dependence accordingly of 



0.842 



1.167 + (r/r s ) 3 



(42) 



This leads to a gradual flattening of the spectral index with 
r. The quantitative behaviour of this simple approximation 
is in agreement with the data shown for the few supernova 
remnants for which we have adequate distances or know their 
individual time evolutions (Anderson et al. 1991, Fig. 1). 

We note that these effects depend on the ambient density 
n such that for a low density environment the Alfven velocity 
is high and hence second order Fermi acceleration may then 
become important. However, we have to emphasize that second 
order Fermi acceleration depends on the average strength of 
the Alfvenic wave field across the entire shell region, and since 
that is unknown, we can derive the temporal change of the 
spectral index only to within an unknown factor, which might 
be rather small. Thus we see two ways to explain the trend for 
overall spectral index to flatten with diameter (Berkhuijsen 
1986, Anderson et al. 1991): 

First, there is the possibility that for larger and thus older 
remnants the spectral index may in fact be due to a combi- 
nation of various components that may have formed over the 
long time scale available: 1) HII regions with their rather flat 
radio spectrum from free-free emission, 2) a pulsar driven su- 
pernova remnant which also has a flatter radio spectrum, and 
3) an underlying synchrotron spectrum of a shell-type super- 
nova remnant as discussed in this paper. Furthermore, the su- 
pershells often found in other galaxies require the energy input 
of many stellar winds and/or supernovae and might be an anal- 
ogy to the extremely large supernova remnants detected in our 
Galaxy, and not just attributable to a single event. 

Second, one could argue that only a few remnants sur- 
vive to the age where second order Fermi acceleration becomes 
important, and so these remnants are not significant for the 
typical injection spectrum for the release of cosmic rays into 
the interstellar medium. The distribution of supernova rem- 
nants is given by a r +1 ' 2 power law in equal logarithmic bins 
of observed radius, for any given environmental density n . 
The adiabatic expansion covers a range of only 8.1 nZ 0,12 from 
free expansion to cooling and break up, and thus this range is 
rather insensitive to the density. However, the breakup radius 
depends as rio'^ 11 on density and so the typical range of at 
least a factor of 100 in the interstellar medium density (see, 
e.g., Fig. 12 in Berkhuijsen 1986) likely to be experienced by 
exploding stars will shift the distribution around by an extra 
factor of ten, totally mixing the diagram (see, e.g., Fig. 11 in 
Berkhuijsen 1986). Thus, it is close to impossible to disentan- 
gle these various effects all of which lead to selection effects in 
any observational survey. We suspect, however, that it is this 
second alternative which is more reasonable in our Galaxy. 



6.1. The time evolution of the radio emission 

Calculating now the radio emission contribution from the shell 
of 7" t/2/(3 C/i), we note that the ram pressure is proportional 
to Ui, which in the Sedov solution [Eq. (2)] is proportional to 
r -3 . The total radio emission from the shell is then constant 



T c A in23 D 1.710 ci -0.710 -Itt -1 

L v = 5.410 S_ 5 .3 77-1 #51 vg. ergs Hz 



(43) 



and so does not depend on the interstellar medium density; 
since the adiabatic expansion is a condition of constant energy, 
all such dependences drop out. We have used here as reference 
0.1 for 77, 5 fiG for the unperturbed interstellar magnetic field 
strength, and 1 GHz for the emission frequency observed. We 
note that the emission here is coming from a shell of thickness 
r/12; additional emission further inside might arise if there 
is a reverse shock resulting from the transition of free expan- 
sion to adiabatic expansion. Evidence for reverse-shock heat- 
ingof ejecta may exist in X-ray images of Tycho (Seward et al. 
198J ), and has also been c onstrue d from X-ray spectroscopy of 



SN 10 06 (Hamilton et al. |l986a|) and Tycho (Hamilton et al. 
1986b[ ). 

Furthermore, there is observational evidence that the syn- 
chrotron emission systematically samples the higher magnetic 
field substructures, by a factor of up to 5 and so we can 
expect the synchrotron emission to be higher by factors up 
to 15.7, leading to an implied radio luminosity of up to 
8.5 10 24 ergs _1 Hz _1 , everything else being equal. This lumi- 
nosity is in agreement with the upper limit of the distribution 
of observed luminosities, suggesting that in their initial phases 
the efficiency of putting kinetic energy into an energetic elec- 
tron population is indeed of order 0.1. 

However, when we derive the normalization factor for the 
particle spectrum to use in the approximation that the total 
energy in all relativistic particles is 77 times the flow energy, we 
have implicitly used a constant particle spectrum. After having 
shown how we might understand the early time evolution of the 
spectral index, we now have to calculate how this influences the 
normalization; since with an increasingly flatter spectrum the 
energy in the particle population shifts to higher values, there 
are increasingly fewer particles at lower energy, the very parti- 
cles responsible for radio emission. However, this effect is only 
strong for relativistic particles and is counterbalanced by the 
non-relativistic sector of momentum space, so that the entire 
effect is only a factor of order 2 for a change from a synchrotron 
spectral index of —0.7 to —0.8. Therefore, the radio emission 
is expected to decrease slowly with time. As soon as electron 
injection ceases, then simple aging will decrease the electron 
spectrum rapidly and Shklovsky's argument holds (Shklovsky 
1968, §7), demonstrating how the electron population decreases 
simply due to overall adiabatic losses. 

7. The injection into the interstellar medium 

Since the particle spectrum flattens with time and the observed 
cosmic ray spectrum corresponds to the early spectral index 
of Sedov expansion explosions, there is an apparent discrep- 
ancy. This discrepancy can be resolved by noting that a Sedov 
type expansion breaks up when the shell becomes unstable to 
cooling; the observations and our interpretation then suggest 
that most supernova remnants break up and release their en- 
ergetic particle population into the interstellar medium long 
before second order Fermi acceleration becomes relevant. The 
Vela supernova remnant may be an example of a shell that 
has been breaking up for some time. Spectroscopy has shown 
that the motions in Vela are very chaotic (Jenkins et al. 1984), 
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with large departures from simple radial motion, which may 
be evidence for this. 

The evolution of a supernova remnant (see, e.g., for the 
basic physical concepts employed Cox 1972) is then divided 
into four phases: 

First, there is free expansion until the interstellar medium 
overrun by the explosion shock is of about the same mass as 
the mass ejected. This happens at a radius of 



Re 



(44) 



where M c j is the ejected stellar mass and n a is again the 
ambient density in cm~ 3 . Numerical simulations demonstrate 
that the expansion is noticeably slowed down long before this 
critical radius is reached. 

Second, we have a Sedov phase with adiabatic expansion 
and steady fresh electron injection. The observations of radio 
emitting stars suggest that there is a critical velocity for a 
shock, below which electron injection is severely restricted. The 
critical speed is of order 



Ui, e 



c(m c /m p ) 1/2 ~ 7000 kms~ 



(4£ 



as suggested by the radio luminosities of OB and WR stars 
on the one hand, and radio supernovae on the other (CR II); 
this comparison suggests just the approximate velocity, not the 
actual number, but the particular velocity given here is consis- 
tent with the stellar data. It is highly unlikely that the same 
critical velocity would be relevant in the interstellar medium, 
and so we only use the concept that there is such a critical ve- 
locity and adopt as reference 1000 km s" 1 (= C/ C rit,-2.5, in units 
of c) on the grounds that this would correspond to roughly the 
same Alfvenic Mach number as the supernova shocks in stellar 
winds as argued in papers CR I and CR II without any physical 
justification at this stage. We plan to test such a concept with 
further stellar examples of nonthermal radio emission. In our 
picture the injection is done by both the bow shocks around 
the fast moving knots and the overall shock. Hence the critical 
radius is given by 



Hcrit = 9-0 (En/no) 1 ' 3 Ujf* 



pc. 



(46) 



Third, when the expansion velocity has decreased below 
the critical speed, the electron injection ceases and we have 
further adiabatic expansion, but with the electron population 
only modified by adiabatic losses. Adiabatic losses of a rela- 
tivistic particle decrease its energy by the ratio of the initial 
and final radius scales. 

And finally, fourth, we have the break-up phase, which oc- 
curs because the cooling layer is thermally unstable (McCray 
et al. 1975, Chevalier and Imamura 1982, Smith 1989) leading 
finally to a break-up of the supernova shell. The cooling radius 
is (Cox 1972) 



Rc 



15.6 E^ 11 (10 22 L)- 2/n 



-5/11 

n pc, 



(47) 



where L is the cooling coefficient; the cooling curve has 
been compared with recent detailed calculations by Schmutzler 
& Tscharnuter (1993) and is still a very good approximation 
at the temperature range of interest here. During this time 
the energy of any individual relativistic electron decreases by 
adiabatic expansion as the ratio of the radii from the time 



when injection ceases to the time of break up. Thus the energy 
density of the electron population (from conservation of the 
adiabatic moment only and disregarding spatial dilution, see 
below) decreases by the factor (p = 2.420) 



' RmolV '= 2.18 J B 5 - 1 09 n o - 17 (10 22 L)- - 26 C/ c r f t B _ 2 . 5 , (48) 



showing a very weak dependence on the environmental den- 
sity n a . The dependence on the critical velocity is very nearly 
linear. The simple spatial dilution of the energy density is the 
same for both protons and electrons and drops out in the ra- 
tio, as long as we are in the Sedov phase (see Shklovsky [1968], 
eq. 7.27). For a tenuous interstellar medium of H-atom number 
density about 0.01 cm , the factor is increased to nearly 4.8. 

This intermediate switch from electron injection with 
steady acceleration to a simple adiabatic loss regime will de- 
termine the net scaling of the power of the electron population 
to that of the proton population in the cosmic rays. The ob- 
servations suggest that from 1 GeV the electron number is 
only about 0.01 (Wiebel 1992) of the protons. From this ob- 
served ratio the energy density ratio integrated over the entire 
relativistic part of the particle spectrum, protons relative to 
electrons, for the injection spectral index of —2.42, is given by 
4.3. Assuming then that the ratio of electrons to protons is not 
influenced by propagation effects, this suggests an expansion 
of a factor in radius of order 3 between the time when electron 
injection ceases and the time when cooling breaks up the shell 
of the supernova remnant; here we assume that electrons and 
protons originally have comparable energy densities of their 
relativistic particle populations. 

Interestingly, this is very close to the ratio implied above 
for normal parameters. Below, where we argue about the max- 
imum particle energy that can be attained, we see that the ten- 
uous part of the interstellar medium is likely to be the phase 
responsible for the injection. This is then fully consistent with 
our adopted particular value for the critical velocity of 1000 
km s _1 and suggests that we may have identified the physical 
reason for the observed electron/proton ratio. 

8. The maximum energy of particles 

The maximum energy particles can possibly reach is given in 
Sect. 3, and depends linearly on the magnetic field. Near the 
symmetry axis we have diffusion parallel to the magnetic field, 
and, given the strong turbulence induced by the shock we take 
there the Bohm limit for the upstream diffusion coefficient: 



1 E 



3 ZeBi 



(49) 



At the equator we use our radial upstream diffusion coeffi- 
cient derived earlier, here for a strong shock, of 



K± = 



rUi 
12 



(50) 



Combining Eqs. (49) and (50) for the two upstream diffu- 
sion coefficients (Jokipii 1987) to obtain an acceleration time 
scale and setting that equal to the upstream flow time scale 
of r/(3Ui) then yields a condition on the maximum particle 
energy that can be reached. This condition is 
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r 



4 E 
3 ZeB 1 

This leads to 



c 



3Uj 



(l-^ 2 ). 



4 c 



(51) 



(52) 



A more detailed consideration of the strength of the mag- 
netic field as a function of latitude does not change this result. 
In addition, since here the actual coordinate n does not appear 
in the expression, it is obvious that this condition holds locally 
over any area of the spherical shock with a locally homoge- 
neous magnetic field. Therefore, the overall inhomogeneity of 
the magnetic field does not change this result. 

Using normal interstellar magnetic fields of 5 /j,G, and a 
standard Sedov explosion with an arbitrary radius then yields 
a maximum particle energy of 

£ max (protons) = 1.0 10 5 (r/pc)" 1/2 (E 5 i/n ) 1/2 GeV (53) 
and 

£ ma x(iron) = 3.0 10 6 (r/pc)~ 1/2 (F 5 i/n ) 1/2 GeV. (54) 

At breakup, when we suggest the actual mixing of the parti- 
cle population into the interstellar medium occurs, these max- 
imum energies are 



.Emax (protons) = 2.610 For GeV, 
and 

£ ma x(iron) = 7.4 10 5 For GeV. 
where we use 



TP ,7,0.364 / 1n 22 n 0. 091 -0.273 

rcR = F 51 (10 L) n a 



(55) 



(56) 



(57) 



The protons ought to dominate due to the normal abun- 
dances in the interstellar medium, and higher energies for pro- 
tons are clearly reached in the tenuous hot part of the inter- 
stellar gas where the densities are of order 0.01 cm~ 3 , and so 
maximum proton energies are possible to about 1.0 10 GeV. 
These numbers are similar to earlier estimates of the maximum 
particle energy in a Sedov expansion phase of a supernova. In 
paper CR IV we obtain from a fit to the air shower data an 
estimate for this number of 1.2 10 5 GeV, fully consistent with 
these arguments, and suggesting that the cosmic ray injection 
is most effective in the tenuous part of the interstellar medium, 
an argument which has been made before on quite different 
grounds. 

9. Summary 

In this third paper of a series on the process of cosmic ray ac- 
celeration we have again used a basic conjecture (Sect. 3) on 
the diffusion of particles in a shock perpendicular to the mag- 
netic field. In a similar vein we have introduced a number of 
heuristic arguments which require testing against observations. 
We have already tested and further explored the consequences 
of this concept in another paper (CR II) on the nonthermal 
radio emission of Wolf-Rayet stars and demonstrate that our 
concept can produce the proper radio spectral indices, lumi- 
nosities and temporal behaviour. In a further communication 



we have tested in detail the predictions of this model as re- 
gards the chemical abundances of cosmic rays (CR IV). What 
we have shown in papers CR I, II, IV and here, is that the entire 
cosmic ray spectrum from 1 GeV out to 3 10 9 GeV, including 
the knee energy and the chemical abundances, can be under- 
stood with the same physical ingredients already well-tested in 
the solar wind shock, namely a strong shock, diffusion, drifts, 
convection, adiabatic cooling, and injection history. 

In this paper in particular, we have used the concept intro- 
duced in paper CR II of a critical velocity for electron injection 
to fashion an argument for the observed electron/proton ratio 
in cosmic rays. We have derived an estimate for the upper par- 
ticle energy that can be reached in such explosions, of order 
10 5 GeV, close to what had been argued previously and what 
is required by the air shower data (paper CR IV) . We have also 
derived the proper cosmic ray spectrum, consistent with direct 
measurements as well as with radio observations. It obviously 
follows that the abundances in such cosmic rays ought to be 
close to those of the interstellar medium itself, as is observed. 
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